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_H ^ Abstract. We present new high-sensitivity high-resolution mm-wave observations of the Onsala 1 ultra-compact 

r^ ■ HII region that bring to light the internal structure of this massive star forming cloud. The 1.2 mm continuum 

map obtained with the TRAM 30-m radiotelescope (~ll" resolution) shows a centrally peaked condensation of 1' 

size (~0.5pc at the assumed distance of 1.8 kpc) which has been further investigated at higher resolution in the 

3 mm continuum and in the emission lines of H^^CO^ J=l-0 and SiO J=2-l with the IRAM Plateau de Bure 

interferometer. The 3 mm data, with a resolution of ~5"x4", displays a unresolved continuum source at the peak 

of the extended 1.2 mm emission and closely associated with the ultra-compact HII region. The H^^CO"'' map 

^ ' traces the central condensation in good agreement with previous NH3 maps of Zheng et al 1985 However, the 

\^ I velocity field of this central condensation, which was previously thought to arise in a rapidly rotating structure, 

^^ . is better explained in terms of the dense and compact component of a bipolar outflow. This interpretation is 

>0 ' confirmed by SiO and CO observations of the full region. In fact, our new SiO data unveils the presence of 

^y , multiple (at least 4) outfiows in the region. In particular, there is an important center of outfiow activity in the 

region at about 1' north of the UCHII region. Indeed the different outfiows are related to different members of 

the Onsala 1 cluster. The data presented here support a scenario in which the phases of massive star formation 

begin much later in the evolution of a cluster and/or UCHII region last for much longer than ICyrs. 

^\ 
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CZJ 1. Introduction et al ll998l Watson et al. I2003|l . We therefore adopt the 

f^ near-distance of 1.8 kpc to this source for the purposes of 

^ _ Onsala 1 (ONI) is an ultra-compact HII (UCHII) region |.^jg work 

• 'H . situated in the densest part of the Onsala molecular cloud 

rS [ (Israel & Wootten ll983(l . It is one of the most compact Several masers such as OH (Ho et al. 11983 ). H2O 

j^ . and isolated UCHII regions known to be associated with (Downes et al. ll979|l . CH3OH (Menten ri991|l and HCHO 

■ " " ' several signposts of massive star formation such as dense (MacLeod et al. 11998(1 are known to coincide with ON 1 . 

gas, high velocity gas, masers, and near-infrared sources. Spatially separated blue and red shifted OH maser spots 

However, a single luminous (1O''.2L0) far-infrared (FIR) around ONI have been identified by Desmurs & Baudry 

point source namely IRAS20081+3132 is associated with |)1998|l using VLBI techniques. Observations using dense 

this UCHII region. The luminosity estimated from FIR gas and dust tracers have shown a compact, nearly circu- 

data and the centimeter (radio) emission agree weU, within lar dense core surrounding ON 1. Maps of 350^m contin- 

a factor of 100 L0, and indicate a spectral type of BO. 3- uum emission (Mueller et al. l2()()()|l and CS J=5-4 emission 

BO. 6 (MacLeod et al. 11998(1 . At a spatial resolution of (Shirley et al. 120031 show dense circular condensations 

0.7"-1.5", the UCHII is detected as an unresolved source around the source. Emission from transitions of CH3CN 

at 2 cm & 3 cm wavelengths (Kurtz et a UBMj) . Higher an- J=12-ll (Pankonin et al. IMITIl and HNCO (Zinchenko 

gular resolution(0.1") observations at 1.2 cm resolve the et al. '2000) are also found to peak on ONI. Using VLA 

UCHII region into an approximate shell structure of size observations at ll"resolution in the (J,K)=(1,1) hue of 

0.6"(Turner & Matthews CnHl) • This source is placed at NH3, Zheng et al. iP^ hereafter ZHRS85) argued for 

a distance of 1.3-6 kpc by various authors in the htera- the presence of a rapidly rotating gas condensation cen- 

ture. However, a distance of less than 2 kpc is favored by tered on ON 1. However, the regions of massive star forma- 

most authors for many valid reasons (see e.g. MacLeod tion often possess an intricate kinematic behavior involv- 

ing motions of rotation, infall, and outflow, and the actual 

Send offprint requests to: M. S. Nanda Kumar existence of a rapid gas condensation in ON 1 needed con- 
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firmation with higher resolution observations. In this pa- 
per, we present high-sensitivity high-resolution (~5") line 
and continuum observations which are used to investigate 
in detail the dense core and outflowing gas associated with 
ONI. 



2. Observational Data 

The IRAM 30 m radiotelescope at Pico Veleta (near 
Granada, Spain) was used to map the A 1.2 mm contin- 
uum emission from the ON 1 region in December 2000. 
We used the MAMBO 1 bolometer array to produce an 
on-the-fly map with a scanning speed of 4" s^^, a wobbler 
period of 0.5 s, and a wobbler throw of 53". The central 
frequency of the detectors was ^^240 GHz, and the band- 
width '^TO GHz. Sky dips made immediately before and 
after the map were used to correct for the atmospheric 
extinction, and a global calibration factor of 15000 counts 
per Jy (based on an observation of Uranus) was applied 
to the data. The resulting map (produced with the IRAM 
NIC software) has an angular resolution of 11". 

High resolution interferometric observations of ON 1 
region were carried out using the IRAM Plateau de 
Bure (PdB) interferometer (near Grenoble, France) in 
its D (5 antennas) and C (6 antennas) configurations. 
Observations using the compact D array configuration 
were carried out on the nights of 28 May 2002 and 
20 July 2002 for two adjacent fields; one field centered 
on the ON 1 FIR point source/UCHII region and an- 
other at offset I'north. The C array configuration ob- 
servations for the same fields were obtained during the 
nights of 9 and 10 January 2003. Each field was a mo- 
saic of 9 pointings arranged in a 3 x 3 matrix with a spac- 
ing of 12" which is roughly equal to half the primary 
beam size. We used the SIS heterodyne receivers both 
at 87 GHz and 230 GHz. The 230 GHz data was mostly 
unusable due to summer weather conditions. Six spec- 
tral correlators were used to select three spectral lines, 
namely, SiO J=2-l(:y=86. 846891 GHz), Hi3C0+ J=l-0 
(zy=86. 754330 GHz) and CO J=2-l (i^=230.538GHz) and 
two continuum bands centered at 1.2 mm and 3 mm wave- 
lengths. The system temperature was typically 100 K at 
3 mm and >400K at 1.2 mm with a phase noise of ~25°- 
35°at 3 mm and >40°-50° at 1.2 mm. The D configura- 
tion provided baselines of 24m-82m while the C config- 
uration resulted in baseline range of 48m-229m. We used 
the sources 2023-1-336 and 2013-1-370 as the primary phase 
calibrators and MWC349, 2013-1-370 and 3C345 for fiux 
calibration. Synthesized clean beams at 3 mm measured 
'^5.35"x3.7" (PA=56°) which represent the typical res- 
olution of the interferometric data presented here. The 
resulting maps does not include any short spacing data 
from single dish observations and thus filter out all large 
scale structures. 
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Fig. 1. 1.2 mm continuum emission (dotted contours) 
mapped with IRAM30 m telescope overlayed on a K-band 
image (Kumar et al. 12003(1 of the Onsala 1 region. Contour 
levels are 1,2,3,6,9,12,15,18 & 21 Jy/beam and has an rms 
of 140mJy/beam. Thick contours enclosing grey shades 
represent the 3 mm continuum emission mapped with the 
IRAM PdB interferometer. Levels are 0.1, 0.15 & 0.2 
Jy/beam (rms = 1.6m Jy/beam). The star symbol rep- 
resents the position of the FIR sources and the UCHII 
region. Solid lines represent blue-shifted and dashed lines 
represent red-shifted H^^CO"*" J=l-0 emission, integrated 
in the velocity range V;sr±4.5kms^^. In both cases, the 
contour levels start at 0.5Kkms^^ (rms = 0.06Kkms~^) 
and increase in steps of 0.5Kkms~^. The beam sizes of 
the 1.2mm data(dotted circle) and the 3mm data (bold 
ellipse) are shown enclosed in grey shades at the bottom 
right corner of the figure. 



3. Results 

3.1. Continuum observations 

In Fig.l we depict the results of observations that repre- 
sent the dense gas and dust in the region of ON 1 . The star 
symbol shows the coinciding position of the UCHII region, 
the luminous IRAS point source and all the maser spots 
that represent the massive young stellar object, which we 
shall henceforth refer to as ONI. The 1.2mm continuum 
emission is extended over a size scale of ^70" and elon- 
gated in the North-East direction. The 1.2 mm continuum 
emission is a sensitive tracer of the warm dust and there- 
fore indicates the dense cloud surrounding ON 1. The spa- 
tial extent and morphology of the 1.2 mm continuum map 
is very similar to that of 350/im map (Mueller et al. l2000|l . 
the CS emission line maps (Shirley et al. I2()03|l and also 
the integrated NH3 emission (ZHRS85) which has been 
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thought to represent a rapidly rotating condensation. The 
source ONI and the surrounding embedded cluster all 
appear well immersed in the observed 1.2 mm emission. 
Therefore the 1.2 mm emission region in Fig. 1, which is 
representative of the condensation unveiled by all other 
dense gas tracers, can be considered as the dense core 
in which the ON 1 cluster has been forming. The 3 mm 
continuum emission reveals an unresolved source centered 
on ONI. The interferometric observations filter out the 
extended emission and show the underlying small scale 
structure. Therefore the 3 mm continuum emission cen- 
tered on ON 1 is indicative of the major condensations 
within the more extended clump traced by the 1.2 mm 
emission. However, 10/xm observations (Kumar et al.[2003 ) 
reveal two closely situated sources coinciding with the un- 
resolved 3 mm source indicating that ON 1 is likely to har- 
bor a binary protostar. 



The 1.2 mm emission maps can be used to estimate 
the mass of the ON 1 clump and the 3 mm continuum 
emission the mass of the unresolved central source(s). The 
black body fits to the Far-Infrared (FIR) flux density dis- 
tribution of ON 1 indicates a dust temperature of 57 K 
(Mueller ct al. I2()()()|l . Since the peak of the dust conden- 
sation coincides with an UCHII region, this value is likely 
an upper limit. If we assume a dust absorption cocfh- 
cient of Ki.2mm=0.0005kg^^m-^ (Henning et al. .1985,1 and 
a distance of 1.8 kpc to the source, the observed 1.2 mm 
continuum emission represents a cloud mass of ~386OM0 
for T=50K and 5000Mq for T=40K. The source aver- 
aged flux density (3.6 Jy/beam) imply hydrogen column 
densities of 1.4xl0^''cm~^ resulting in Ay ^^700 mag. The 
peak column density is even higher at 8.2xl0^^cm~^ im- 
plying a few thousand magnitudes of visual extinction. 
However, these numbers are to be treated with caution 
due to the presence of the UCHII region. The observed 
emission does not faithfully represent the temperature of 
the warm dust since the heating effects of the UCHII 
region can contribute to a majority of emission close 
to the peak. The 3 mm emission from the unresolved 
source represents a mass of ~400Mq, where we have used 
'«3TOm=6.6xlO~^kg~^m^. This value of K^mm is obtained 
by assuming a spectral emissivity index (3—2. There is 
increasing observational evidence which suggest that the 
value of (3 is likely much lower (/? '^l) around massive 
protostellar sources (ex: Hunter et al l2()()()|l . If we assume 
/? ~1 then the resulting estimates of the embedded mass 
in the observed 3 mm emission will be lower than quoted 
above. In any case, the derived mass from 3 mm emission 
is probably an over estimate since 3 mm continuum emis- 
sion is known to trace free-free emission also and may not 
be a good tracer of the dust alone (Downes et al. ll992|l . In 
view of these facts we can safely assume that the enclosed 
mass in the unresolved central object is not greater than 
--3OO-4OOM0. Thus about 6-8% of the overall dense core 
mass seems to be in the form of protostellar material at 
the center. 



3.2. H^^CO+ J=l-0 Emission 

In Fig. 1 solid and dashed contours show blue and red 
shifted H^'^CO"'" J=l-0 emission, respectively, integrated 
in the velocity range V;sr±4.5kms~^. As evident from 
the figure, the H-'^'^CO"'" emission is centered on the source 
ON 1 and reveals a spatially resolved bipolar structure. 
The position angle (PA) of the axis of this well-collimated 
bipolar structure is ~44° and roughly coincides with 
the elongation displayed by the 1.2 mm contours. The 
H^^CO"*" emission also coincides well with the NH3 con- 
densation mapped by ZHRS85. The projected length of 
this bipolar emission is ^--^30" corresponding to about 
0.27 pc at the assumed distance of 1.8 kpc. Each lobe of 
this compact feature appears to consist of two distinct 
clumps. 

In Fig. 2 we show a velocity-position map obtained 
along the axis (PA~44°) of this bipolar feature which re- 
veals a structure that is distinctly different from earlier 
claims of ZHRS85. The vertical line represents the sys- 
temic LSR velocity at 12kms~^. The red-shifted south- 
ern lobe and the blue-shifted northern lobe can be seen 
clearly separated into two components on the PV dia- 
gram. A main component of the emission (denoted by A 
& B in Fig. 2) stretches to a radius of '^5" and a satel- 
lite component comprising of two knots (denoted by C & 
D) extends to a radius of ~10". These 4 clumps corre- 
spond to the clumps found on the images (Fig.l). As can 
be seen from Fig. 2, all of the found emission is spread 
along a straight line on the PV-plane showing increasing 
velocity with increasing distance from the central position 
following a Hubble law. We have not subtracted the bright 
continuum emission from this data. The continuum emis- 
sion can therefore be seen close to the source (offset = 0) 
extending on cither side of the Yisr, up to the limits of 
our observing setup with the spectrometers. 

Since the H^'^CO'^ J=1~0 emission is symmetrical both 
in position and velocity with respect to the central posi- 
tion and velocity, it is natural to conclude that such emis- 
sion arises in a bipolar outflow emerging from the central 
object. Note that the observed kinematical structure can 
not be easily explained with a rotating disk. First of all, 
the total size of this structure (>0.2 pc) is too large to be a 
disk and, moreover, a disk would produce just two peaks of 
emission in the PV diagram, not four peaks as the A,B,C, 
& D peaks which are actually observed in Fig. 2. The 
Hubble-like velocity law distribution which is observed 
in Fig. 2 (velocity proportional to the distance) is also 
very typical of young bipolar outflows. A straight line fit 
through all the knots shows a gradient of '^30kms~^ pc~^ 
at the assumed distance of 1.8 kpc. These velocity limits 
and the velocity gradient are similar to the values obtained 
from NH3 observations by ZHRS85 except that the gra- 
dient of ^llkms~^ pc""'^ estimated by ZHRS85 is due to 
their adopted distance of 3.5 kpc. The interpretation that 
the H^'^CO^ emission arises in the compact component 
of a bipolar outflow is confirmed by the SiO high- velocity 
data shown in the following section. 
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Fig. 2. Velocity-Position diagram for the H^^CO+ J=l-0 
emission shown in Fig. 1, obtained along the axis of the 
bipolar emission at PA~44°. The contours levels start at 
1 K and increase in steps of 1 K. 



Assuming standard LTE conditions, the observed 
H^'^CO"'" emission can be used to compute the column den- 
sities and enclosed mass. If we assume a dipole moment of 
Ai=4.08debye (Haese & WoodflnZ^and an Hi^CO+ZHs 
ratio of 3.3x10"" (Blake et al. HMTl in OMCl), the 
observed fluxes indicate an enclosed mass of ~28OM0- 
33OM0for temperatures of 40K-50K. The estimated mass 
can be uncertain by a few hundred solar masses depending 
on the assumed value of H^'^C0+/H2 ratio and excitation 
temperature. In any case, the estimated H^^CO"'" mass is 
of the same order than the mass evaluated from the 3 mm 
continuum emission. 



3.3. SiO & CO observations 

In Fig. 3 we present the results of the SiO J=2-l and 
CO J=2-l high-resolution observations in all the region 
around ON 1 including the north zone in which the 1.2 mm 
emission is significant. The SiO emission is found at ve- 
locities of V;sr±4-5kms~^, with respect to the systemic 
velocity, whereas the CO emission appears to extend up 
to Vjsril2kms~^. The CO J=2-l data at 1.2mm is quite 
poor due to summer weather conditions, but the knots 
shown in Fig. 3 have high signal to noise ratio (10-15(t). 
Based on complementary blue and red shifted components 
we can identify at least four main outflows in the region. 
However, we caution that these maps do not reveal most 
of the large scale structures in the region, since we do not 
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Fig. 3. SiO J=2-l (thin lines) and CO J=2-l (thick hues) 
integrated emission overlayed on a continuum subtracted 
H2 narrow band image at 2.122/.tm. SiO emission was in- 
tegrated with a velocity interval of V/sr±4kms^^ and 
CO emission was integrated in the range V/sr±12kms~^ 
Solid lines represent blue-shifted and dotted lines rep- 
resent red-shifted emission. SiO contour levels start at 
0.2Kkms^^ (rms = 0.06Kkms~^) and increase in steps of 
0.15Kknis^^. CO contours levels start at 2Kkms~^ (rms 
— 0.22 Kkms) and increase in steps of IK km s~^. The star 
symbol represents the position of the FIR sources and the 
UCHII region as in Fig. 1. The dotted straight line repre- 
sents the axis of the H^^CO^' outflow. The beam sizes of 
the 3mm and 1mm data are shown by filled ellipses at the 
bottom right corner of the figure. 



have any short spacing data from single dish telescopes 
included in the maps. 

The first outflow is seen close to the central FIR source. 
The most internal component of this outflow is at PA '^45° 
in good agreement with the H^^CO+ outflow component 
discussed in Section 3.2. As for the H^^CO"*" data, the 
blue-shifted emission is placed at the NE of the source, 
while the red-shifted emission is placed at the SW. We 
then beheve that the SiO and the H^^CO+ high veloc- 
ity emission is tracing the same compact bipolar outflow 
which is emerging from the very vicinity of the FIR object. 

Most remarkably, the region placed 40" north from 
ON 1 FIR appears to be an important center of outflow 
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activity: a prominent high-velocity SiO feature which is 
marked as Atiue in Fig. 3 stretches approximately 30" 
in the north-south direction forming a kind of highly- 
coUimated jet. the red shifted counterpart to this jet-like 
structure is not so obvious. We think that the bright red 
knot (Ared) on the south can be associated with Atiue be- 
cause both structures are well aligned and because there 
is no other equally prominent blue component in the re- 
gion. If these two features Ared and Auue were actually 
part of the same outflow, this would have a total size of 
at least 1' or 0.5 pc. We caution however that it is difh- 
cult to identify bipolar counterparts at large scale in such 
complex regions. 

A smaller outflow denoted by B shows both blue and 
red shifted lobes and extends to about 25" in the east-west 
direction. This flow is traced even by the poor quality 
CO data appearing as spatially separated blue and red 
knots. The driving source of this outflow has not been 
identified neither in the near-infrared nor in the mm-wave 
observations. The 3a noise limit of 5 mJy of our 3 mm 
continuum maps imply a limiting detectable mass of 5- 
6 M0 . Therefore this outflow is likely driven by a low mass 
object. 

Finally, a short east-west jet-like feature, denoted by 
C, is emerging from the dense blob of SiO emission sur- 
rounding the ON 1 source represents another distinct bipo- 
lar outflow. Indeed, this flow is traced more prominently 
by the symmetrically placed CO knots in an axis that 
closely matches with the axis of the overall bipolar distri- 
bution of near-infrared H2 emission. It is unclear if this 
outflow arises from the same object than the dense and 
compact NE-SW outflow seen in H^^CO"*" and SiO, and 
which has been discussed above. Since the FIR object is 
probably a binary system (Kumar et al. 12003(1 . it may 
be not surprising to find two different outflows emanating 
from its vicinity. 

It appears anyway that high velocity emission is spread 
all over the cluster forming molecular cloud. Although our 
observations indicate the presence of real multiple out- 
flows, it is also likely that some of these outflows have 
interacted with the surrounding cluster and broken into 
multiple components. This may be particularly true for 
the east-west flow which is traced by the CO knots, a 
faint SiO jet and largely by the NIR II2 emission. 

4. Discussion 

ON 1 was described to be a rapidly rotating condensation 
by ZHRS85 based on NII3 observations. The higher reso- 
lution data presented here along with other dense gas and 
dust tracing observations indicate that the NH3 conden- 
sation described by ZIIRS85 actually corresponds to an 
outflow. The H^'^CO^ observations discussed in Section 
3.2 has the orientation and velocity ranges similar to the 
NII3 features. However, for the first time, our observa- 
tions shows a spatially resolved bipolar feature whose PV 
diagram is typical of outflows. Perhaps the H^'^CO^ ob- 
servations alone would not have resolved the issue if this 



is an outflow or a rapidly rotating envelope, particularly 
because of the velocity structure close to the source and 
the low values of Aw — V/sr±2-3kms~^. But the SiO emis- 
sion which traces other outflows in the region show only 
slightly higher values of Av = V(sr±4-5kms~^. The sym- 
metrical lobes of H^^CO^ emission indicate that the out- 
flow has a small inclination angle with respect to the plane 
of the sky. This implies small line-of-sight velocities of the 
H^^CO+ structure in comparison to the SiO data. Indeed 
all other dense gas tracers from this region such as CS 
(Shirley et al lM75|l . CH3CN (Pankonin et alEUHII) and 
HNCO (Zinchenko et al lMTOIl display similar V;^^ ~11- 
12kms~^and Av '^4-5kms~"'^. Therefore it is clear that 
the H^'^CO"'" presented here spatially resolves the NH3 
data from ZIIRS85 and reveal an outflow well placed in the 
plane of the sky and possibly slow moving. A similar exam- 
ple where a large scale coherent NH3 structure (Jackson 
et al. I1988|l is aligned well with an outflow (Bachiller & 
Cernicharo nMljl is found in the case of NGC6334I. 

In the light of above results, the following scenario ap- 
pears most plausible for the ON 1 region. The UCHII re- 
gion ON 1 is a relatively isolated massive star forming re- 
gion associated with a small embedded cluster, all of which 
is immersed in a molecular clump of few thousand solar 
masses. The vicinity of the UCHII region appears to drive 
at least two (the H^'^CO'*" compact outflow, the flow traced 
by CO knots and H2) and likely three (SiO flow A) major 
outflows indicating multiple embedded sources of similar 
evolutionary states. While the exact driving sources for 
each of these flows are not clear, it is very likely that the 
H^^CO"'' flow originate from one of the two mid- infrared 
(10/im) sources (Kumar et al. I20()3|l centered on ONI. 
These two mid-IR sources are neither visible at 3.8/im nor 
distinguishable by the millimeter observations. However, 
given that both sources show lO/^m emission and are asso- 
ciated with an unresolved mm continuum peak, and likely 
drive outflows, they may be at a similar evolutionary state. 
Thus the ON 1 source could harbor a massive binary pro- 
tostar. 

The H^^CO^ and SiO outflows display very low line- 
of sight velocities. Low values of line of sight velocities 
in outflows can be expected if they lie in the plane of 
the sky as appears to be case for H^'^CO"'" flow. But the 
SiO flow does not obviously appear to be oriented in the 
plane of the sky therefore indicating that the low line of 
sight velocities are representative of the region. The out- 
flow traced by H^'^CO+ does not seem to show any ob- 
vious SiO counterparts. SiO production requires violent 
shocks and dust grain destruction while H^'^CO+ outflow 
emission only requires the acceleration of the dense core 
material (commonly traced by H-'^'^CO+). Given the high 
density of the central regions of this symmetric clump, 
the outflow traced by H^'^CO^ is likely encountering vis- 
cous forces thus only able to accelerate the gas rather than 
produce violent shocks. Further, since all dense gas trac- 
ers also display similar range of low velocities such as the 
observed outflows, it is suggestive that the outflows are 
experiencing an opposing force in its passage through the 
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very dense, cluster forming core. Since the jet/outflow ex- 
periences the viscous forces uniformly from the point of 
ejection, the force only makes the flow have much lower 
velocities in comparision to other outflows. 

The presence of a surrounding embedded cluster to 
ON 1 UCHII region represents a low mass stellar popu- 
lation that could be anywhere from Class I to pre-main 
sequence phases. If we consider that low mass stars form 
in ~10^yrs and massive stars form in ~10^yrs, the simul- 
taneous presence of two massive outflows and the central 
location of the luminous sources provide two alternative 
scenarios. Either the massive star formation began much 
later than low mass star formation in this cluster form- 
ing core or that the formation of massive star/UCHII 
region phases last for much longer than ~10^yrs. While 
Herbig I|19(i2|l suggested that massive stars form prefer- 
entially late in the evolution of a cluster. Stabler (1985) 
argued that it is indeed so due to statistical reasons. In 
a related context, de Pree et al. ijlH95|l suggest that the 
life times of the UCHII regions can be much longer than 
a few thousand years if they form and exist in dense and 
warm environments. The black body fits to the FIR flux 
densities of several relatively isolated massive star forming 
regions (ex: Hunter et alEnUU! Mueller et al 2000 Beuther 
et al. I2()()2|l suggest high temperatures between 50-lOOK 
and densities >10^cm~'^ supporting the scenario predicted 
by de Pree et al. ( 1995). This is also true for ON 1 in view 
of the high estimated column densities and temperatures. 
However these observations can not effectively rule out 
the possibility that massive star formation initiated much 
later than the surrounding cluster of low mass stars. 

The FIR and radio flux of the ON 1 source and 
the mass estimates from this work suggest that ON 1 
is among the early B type candidate massive pro- 
tostars. This source can be compared to other well 
studied massive protostars such as IRAS05358+3843 
(Beut her et al. IMI^ . IRAS 18556+0 136 (Fuller et al. 
n(M\ or W75N (Shepherd et al. lM?^ . While W75N and 
IRAS18566-t-0136 have UCHII regions, IRAS05358-)-3843 
has none and is likely about the same age as the mas- 
sive protostars such as IRAS20126-I-4104 (Cesaroni et al. 
[TMHll and IRAS23385-I-6053 (Fontani et all2004). Among 
all these, ON 1 is a neat example of a circularly symmetric 
molecular condensation with a centrally placed massive 
binary protostar which is driving multiple outflows and 
surrounded by a young stellar cluster. 
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